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Abstract 

We investigate the static potential in the confinement phase of the SU(2) Higgs 
model on the lattice, where this model is expected to have properties similar to QCD. We 
observe that Wilson loops are inadequate to determine the potential at large distances, 
where the formation of two color-neutral mesons is expected. Introducing smeared fields 
and a suitable matrix correlation function, we are able to overcome this difficulty. We 
observe string breaking at a distance ~ 1.8ro, where the length scale tq has a value 
tq ~ 0.5 fm in QCD. The method presented here may lead the way towards a treatment 
of string breaking in QCD. 
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1 Introduction 



Since the seminal work by Creutz [Q, there have been a number of detailed studies of 
the static potential in non-Abelian gauge theories. At large distances, there is a linear 
confinement potential between a source anti-source pair in the fundamental represen- 
tation of the gauge group. This was clearly established by Monte Carlo calculations of 
the lattice theories close to the continuum limit, both for gauge groups SU(2) and 
SU(3) When these gauge theories are coupled to matter fields in the fundamental 

representation, one does expect that the potential flattens at large distances and asymp- 
totically turns into a Yukawa form. At such distances the potential is better interpreted 
as the potential between static color-neutral "mesons", which are bound states of a 
static color source and the dynamical matter field. So far, this expectation could not 
be verified by Monte Carlo simulations. In particular, in recent attempts in QCD with 
two flavors of dynamical quarks this string breaking effect was not visible The 
same behavior of the potential is, of course, expected in the non-Abelian Higgs model 
in the confinement phase. While the potential in the large distance range could not be 
calculated in early simulations with gauge group SU(2), they yielded some qualitative 
evidence for screening of the potential [lO|Jll|1 . Similarly, the potential between static 
adjoint charges is screened by the gauge fields themselves. Numerical evidence for this 
has been found for gauge group SU(2) [p^| . 

In this letter we consider the SU(2) Higgs model (in four dimensions) as a first 
test case and demonstrate that string breaking exists. In order to roughly compare 
the physical situation in the Higgs model with the one in QCD, we choose a common 
reference energy scale. One immediately thinks of using the string tension. However, 
due to the phenomenon of string breaking itself, the string tension does not have a 
precise meaning (there is no range of r where the potential is linear) and exists only in 
an approximate sense. It is better to fix the overall scale by ro, defined as 



2 F(r) 



= 1.65, 

r=ro 



where F(r) denotes the force derived from the static potential discussed above.Q The 
number on the r.h.s. of eq. (jl|) has been chosen such that ro has a value (estimated 
from phenomenological potential models describing charm- and bottom-onia) of 

r «0.5fm (2) 

in QCD. In the following, all dimensionful quantities are measured in units of ro, which 
is not related to any phenomenological value in the SU(2) Higgs model. In QCD with 
light quarks, the distance r b around which the potential should start flattening off, could 
be estimated in the quenched approximation E 



r b «2.7r (in QCD) . (3) 



1 To be precise: in the theory with matter fields, r 2 F(r) is not monotonic and we expect that there 
are two solutions for ro. The smaller one is to be selected. 
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2 The Higgs model on the lattice 



Working on a four-dimensional hyper-cubic lattice with spacing a we denote by <J>(x) 
a complex Higgs field in the fundamental representation of the gauge group SU(2) and 
U(x,fi) £ SU(2) is the gauge field link connecting x + afi with x. We use the Greek 
symbols, fj,, v, to denote all directions 0,1,2,3 and Roman symbols such as k to denote 
the space-like directions 1,2,3. The Euclidean action for the SU(2) Higgs Model is then 
given as 

S = ^a 4 {^(V M $(x))tV At $(x)-m^t( x )$( x ) + Ao [$t( :r )$( :r )]2} 

X fj, 

+ i t r£/ p }, (4) 

v 

U p {x) = U(x, fj,)U(x + afb, v)U'(x + au, /i)U'(x, v) , = , 

P x 0</i<i/<3 

where we used the lattice covariant derivative, V M $(x) = - [U(x, [j,)<f>(x + afi) — §(x)\. 



Following conventional notation [14], the bare parameters chosen in the present inves- 
tigation are given by 

f3 = 2.2, 

K = (1 - 2A)/(8 - a 2 ml) = 0.274, (5) 
A = k 2 \ = 0.5. 

This point in parameter space is in the confinement phase, fairly close to the phase 
transition, where the model has properties similar to QCD |l4]JTo[| . The lattice resolution 
is of roughly the same size as the one used in the QCD-studies quoted above: from the 
potential calculation detailed in the following section, we obtained 

r /a = 2.83 ±0.03. (6) 

While one cannot expect a very precise result for the detailed form of the potential with 
such a resolution, it is expected that qualitative features are correctly described. 

The string breaking distance depends directly on the mass of the static mesons 
which in turn is sensitive to the bare mass, mo, of the scalar fields. For our chosen value 
of mo, the string breaking distance, r^, is somewhat smaller than the estimate quoted 
above for quenched QCD. Nevertheless, the physical situation is quite similar. We note 
in passing that the self-coupling of the Higgs field, Ao, appears to have little influence 
on the physics of the Higgs model in the confinement phase [|T|]. 

3 Calculation of the potential at all distances 

We now introduce a method, which - as we will demonstrate in the following section 
- allows to compute the potential, V(r), at all relevant distances in the theory with 
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Figure l: The correlation functions used to determine the static potential, from left to 
right: Cwwi Cwm arid Cmm- The lines represent the Wilson lines, the filled circles the 
Higgs field. 



matter fields. Before explaining the details, we would like to mention the basic point, 
which has first been noted by C. Michael [112]. Mathematically, the method is based on 



the existence of the transfer matrix [15] and the fact that it can be employed also when 
external static sources are present (see e.g. [pp. In the path integral, a static source 
at position x, together with an anti-source at position x r = x + rk, are represented by 
straight time-like Wilson lines fixed at these space-positions. These Wilson lines have 
to be present in any (matrix) correlation function from which one wants to compute the 
potential energy of these charges. The space-like parts of the correlation functions, which 
are again Wilson lines when one considers standard Wilson loops, do not determine 
which states appear in the spectral representation of the correlation functions. They 
do, however, influence the weight with which different states contribute. For these 
space-like parts, we therefore use both Wilson lines which will have large overlap with 
string-like states and Higgs fields with a dominant overlap with meson-type states. 
Combining them in a matrix correlation function, the correct linear combination which 
gives the ground state in the presence of charges can be found systematically by the 
variational method described below. 

Let us now give precise definitions of the correlation functions, which are illustrated 
in Fig. |l[ For small values of r or in the pure gauge theory, the static potential can be 
efficiently computed by means of Wilson loops, 



Cww(r, t) 



( tr [W(x, x r ) W(x r ,x r + t6) W^(x + tO, x r + tO) W*( 



X, X 



tO)}) 



t— >oo 



const, x e 



-V(r)t 



(7) 



where x r = x + rk and W(x, y) denotes the product of gauge links along the straight 
line connecting y with x. In addition, one may compute the mass, fj,, of a static meson 
from the correlation function 



C M (t) = ($t(x) W{x, x + tO) $(x + tO)) ~ const, x e _M * . 

f — »oo 



(8) 



Consequently, we expect that for distances significantly larger than ry,, where the rele- 
vant states correspond to weakly interacting mesons, the potential is close to the value 
lim^oo V(r) = 2fj, and can be extracted from the correlation function 

C M M(r, t) = ($t(a; + tO)W^(x, x + i6)$(s) &(x r )W(x r , x r + t6)$>(x r + tO)} . (9) 



3 



In order to investigate all (and in particular the intermediate) distances, we introduce 
a (real) symmetric matrix correlation function Cij(r,t), i,j £ {W, M} with 

Cwm(t, t) = (&(x + t6) W\x, x + t6) W(x, x r ) W(x r , x r + t6) ${x r + t6)) . (10) 



For fixed r, V(r) is extracted from the matrix correlation in the following way [17]. One 
first solves the generalized eigenvalue problem, 

C(t)v a (t,t ) = X a (t,t )C(t )v a (t,t ) X a >X a+1 . (11) 

According to a general lemma proven in [|17]], the ground state energy V(r) = Vo(r) and 
the excited states energies Vx(r), ... are then given by 

aV a (r) = ln(A Q (t- a ,to)/A a (t,t )) + o(e-^W-^( r »*) . (12) 

Here, N is the rank of the matrix C. So far we have N = 2 and for this small value of 
TV the method would require large values of t for the correction terms in eq. (12) to be 
negligible. 

To improve on this, we introduce smeared space- like links [18], setting the smearing 



strength e of [18] to the numerical value e = 1/4. For the Higgs field we employ a 



smearing operator, S, defined as 

5$(x) = P{P*(x)+P U(x,y)^(y)+V £ U(x,y)$(y)} , (13) 



\x — y\ = \/2a \x — y\ — \/3a 



where V§ = <J?/V and U(x, y) represents the average over the shortest link connec- 
tions between y and x. ^ 

For different smearing levels m = l,...,iV$ we evaluate the correlation functions 
introduced above with <&(x) replaced by ^ m \x) = S m <&(x) and similarly with m = 
0, 1, ...,Njj smearing iterations of the space-like gauge fields, where m = corresponds 
to the unsmeared gauge links. This defines a correlation matrix C with rank N = 
+ N v + 1. 

4 String breaking 

We now turn to discuss our numerical results obtained on a 20 4 lattice with periodic 
boundary conditions. The Higgs model can be simulated efficiently with a hybrid over- 
relaxation algorithm. In particular, we implemented the over-relaxation for the Higgs 
field as described in [19|. After a rough tuning of the mixture of the various parts of 



the updating we found that autocorrelations are no problem in our simulations. Errors 
were computed by a jacknife analysis. 



2 The particular form of S was found to be efficient in a study of eq. with various different types 
of smearing operators. Details of this study will be published separately. 
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Figure 2: Ground state energy of a static meson. 



We computed all correlation functions up to r max = i max = 8a ~ 3ro on 2000 
field configurations. We set N$ = 4 and Ny = 2. The variance of the correlation 
functions was reduced in different ways. We used translational invariance to average 
over the base point denoted by x in eqs.([7[jl0|) and cubic symmetry to average over the 
different orientations k. Finally, statistical errors were further reduced by replacing - 
wherever possible - the time-like links by the 1-link integral pQ] , which can be evaluated 
analytically for the gauge group SU(2). 

The first quantity one wants to know is the mass of the static mesons, since this 
fixes the asymptotic value of the potential. It is best computed from the correlation 
Cm extended to a N$ x N$ correlation matrix by considering the smeared 3> fields. As 
for all other energy levels discussed below, the mass is computed from the eigenvalues 
of the generalized eigenvalue problem described above, setting to = 0. In Fig. |2| the 
convergence for large t is exemplified. One observes that the ground state energy can be 
extracted with confidence and with very good statistical precision, a/j, = 0.698 ± 0.002 
(read off at t = 7a which agrees fully with t = 6a). From the correlation function 
without smearing a determination of fi was not possible. 

At all distances, the potential V(r) was then computed using the full 7x7 cor- 
relation matrix. The convergence of eq. ( |l^ ) is shown in Fig. |3| (circles). In earlier 
calculations in QCD, string breaking was searched for in (smeared) Wilson loops. We 
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Figure 3: Comparison of the static potential computed from eq. ( |12| ) using the full cor- 
relation matrix and the sub-block with the (smeared) Wilson loops. Two representative 
values of r are shown. 



studied whether one can succeed in this way by restricting the correlation matrix to 
the corresponding (Ny + 1) x (iVrj + 1) sub-block. The resulting potential estimates 
(triangles in Fig. |3|) are very good at short distances but have large correction terms at 
long distances. Without a very careful analysis one might extract a potential which is 
too high at large distances, when one uses Wilson loops alone. This might explain why 
string breaking has not been seen in QCD, yet. In contrast, from our full correlation 
matrix we can extract safe estimates for V(r) using eq. ( |l2| ) for t ~ 6a. 

We then followed the steps described in || to determine vq and computed the 
potential in units of rp. Considering in particular the combination V(r) — 2/i, one has a 
quantity free of divergent self energy contributions. It is shown in Fig. ^|. The expected 
string breaking is clearly observed for distances r > ~ 1.8 Tq. Around r ~ rt,, the 
potential changes rapidly from an almost linear rise to an almost constant behavior. To 
resolve this transition region using a smaller lattice spacing is an interesting challenge 
to be addressed in the future. 

Overlaps of variationally determined wave functions Vq are a certain measure for 
the efficiency of a basis of operators used to construct the correlation functions. To give 
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Figure 4: The static potential in units of tq. The asymptotic value 2/i has been sub- 
tracted to obtain a quantity free of self energy contributions which diverge like - . 



a precise definition of the overlap, we define the projected correlation function 



Q(t) = VQC(t)v = ^w„e 



-V n (r)t 



with normalization f2(0) = 1 . 



(14) 



Here, n labels the states in the sector of the Hilbert space with 2 static charges. The 



form given above follows directly in the transfer matrix formalism [15]. The positive 
coefficients uj n may be interpreted as the overlap of the true eigenstates of the Hamil- 
tonian with the approximate ground state characterized by vq. "The overlap " is an 
abbreviation commonly used to denote the ground state overlap, ujq. 

We determine loq straightforwardly from the correlation function Q(t) at large t. 
As shown by the circles in Fig. || our operator basis is big (and good) enough such 
that tdQ exceeds about 60% for all distances. It is now interesting to consider also the 
overlap for the Wilson loops alone, i.e. we again restrict the correlation matrix to the 
(JVtj + 1) x (N\j + 1) sub-block. Let us denote the corresponding projected correlation 
function by fi^y an d the overlap by Uq . The computation of loq^ is more difficult, 
because it turns out to be very small at large r. Nevertheless, the expression 



(15) 




Figure 5: Ground state overlaps, once from the full correlation and once from the Wilson 
loops alone. 

converges reasonably fast and cuff can be estimated from the r.h.s. for t ~ 6a. The 
results plotted in Fig. |H| (triangles) show that Wilson loops alone have an overlap which 
drops at intermediate distances and are clearly inadequate to extract the ground state 
at large r. Instead, the full matrix correlation function has to be considered if one wants 
to compute V(r) at all distances. 

5 Outlook 

We have introduced a method to compute the static potential at all relevant distances 
in gauge theories with scalar matter fields. We demonstrated that it can be applied 
successfully in the SU(2) Higgs model with parameters chosen to resemble the situation 
in QCD. There is little doubt that it can be applied for different values of parameters 
in the Higgs model, at least as long as r^/a is not too large. It is then interesting to 
follow a line of constant physics towards smaller lattice spacings in order to check for 
cutoff effects and to be able to resolve the interesting transition region in the potential. 
Furthermore, one might be interested in increasing the mass of the Higgs field in order 
to reach a situation with a larger value of / tq which presumably is even closer to the 
physics situation in QCD. 

From the matrix correlation function one can also determine excited state energies. 
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We have done this successfully for the single meson states but a precise determination 
of the excited potential at all distances needs more statistics. One expects that the 
transition region of the potential can be described phenomenologically by a level crossing 
(as function of r) of the "two meson state" and the "string state" fl2l}| . We are planning 
to investigate this in more detail. So far, we can only say that for r ~ the two levels 
V\{r) and V(r) are close. 

Of course, it is of considerable interest to apply this method to QCD with dynamical 
fermions. It is difficult to predict how well this can be done. Finding good smearing 
operators should not be a problem. However, the correlation functions will involve 
the quark fields and one cannot easily take advantage of translational invariance to 
reduce the statistical errors as was done here. Thus, larger statistical uncertainties are 
expected. On the other hand, in QCD the quark fields are integrated out analytically, 
which usually results in correlation functions with relatively small statistical errors. Also 



new methods |22j should be tried and the final statistical errors can only be determined 
by explicit computations. In any case, the only possible difficulty is expected to be one 
of statistical accuracy. The proper correlation functions can be constructed along the 



lines of ref. Q and of Sect. y. 

The day of completion of this manuscript, a study of string breaking in the three- 
dimensional SU(2) Higgs model appeared |23|]. Both the method applied and the con- 
clusions are very similar to what we find in four dimensions. 

Acknowledgement. We thank Jochen Heitger for helpful discussions. Our sim- 
ulations were performed on the SP2 of DESY at Zeuthen. We thank the staff of the 
DESY computer center for their support. 
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